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ABSTRACT 


Transport  imaging  is  a  SEM-based  technique  used  to 
directly  image  the  motion  and  recombination  of  charge  in 
luminescent  semiconductors^  allowing  for  the  extraction  of 
transport  parameters  critical  to  device  operation.  In  this 
thesis^  transport  imaging  for  ID  structures  was  initiated 
with  work  on  sample  preparation^  modeling  and  initial 
characterization.  One  dimensional  structures  are  being 
integrated  into  forefront  electronics  due  to  their  inherent 
advantages  in  size^  packing  density  and  power  consumption. 

In  this  work  the  one  dimensional  equation  for  steady  state 
minority  carrier  recombination  distribution  solved  for  the 
Gaussian  source  is  derived  and  results  from  numerical 
simulations  are  presented.  The  diameter  of  the  SEM  beam  is 
determined  experimentally  allowing  for  accurate  simulation 
parameters.  Intensity  and  drift  measurements  on  four 
batches  of  top-down  wire  structure  samples^  fabricated  on  a 
AlGaAs/GaAs/AlGaAs  double  heterostructure  using  a  FIB^  are 
presented.  Significant  decreases  in  luminescence  in  FIB 
exposed  regions  are  reported.  Spatial  luminescence  from 
single  bottom-up  GaN  and  ZnO  nanowires  deposited  by  metal 
initiated  metal-organic  CVD  on  Au  and  Si02  substrates  is 
imaged.  CL  spectra  for  GaN  and  ZnO^  with  peak  intensities 
at  3.27  and  3.29  eV^  are  characterized.  Finally^  several 
suggestions  for  further  research  are  offered  including 
transport  imaging  on  contacted  bottom-up  nanowires  and  a 
potential  application  of  transport  imaging  to  FIB  damage 
characterization . 
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I. 


INTRODUCTION 


A.  BACKGROUND 

The  transport  of  free  charge  in  semiconductor 
structures  is  fundamental  to  the  operation  of  many  modern 
optoelectronic  and  semiconductor  devices.  Photodiodes 
including  solar  cells,  photomultipliers  used  in  night 
vision  systems  and  spontaneous  and  stimulated  emission  in 
Light  Emitting  Diodes  (LEDs)  and  lasers  are  just  a  few  of 
devices  and  physical  processes  in  which  free  charge 
transport  parameters  such  as  carrier  mobility  are  of 
critical  importance. 

Since  the  early  1970' s  much  research  in  optoelectronic 
and  semiconductor  devices  has  been  focused  on  decreasing 
device  sizes  due  to  both  practical  advantages,  including 
increased  packing  densities  and  reductions  in  device 
production  costs  and  power  consumption,  as  well  as  physical 
advantages  such  as  dimensional  differences  in  the  density 
of  states  (DOS) [1] .  Quantum  wells,  fundamentally  two 
dimensional  (2D)  structures,  have  found  myriad  applications 
in  the  fabrication  of  laser  diodes,  infrared  (IR)  quantum 
cascade  lasers  and  High  Electron  Mobility  Transistors 
(HEMTs)  to  name  just  a  few  [2]  .  Research  in  the  field  of 
one  dimensional  (ID)  structures,  referred  to  as  nanowires, 
has  already  yielded  many  advances  including  nanowire  field 
effect  transistors  (FETs) ,  crossed  nanowire  p-n  junctions, 
InP  nanowire  LEDs  and  even  ultraviolet  (UV)  lasers  made  of 
single  GaN  nanowires  much  like  the  one  pictured  in  Figure  1 
[3,4,5] .  Accurate  measurement  of  transport  parameters  is  an 
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important  part  of  nanowire  applications  for  all  the  above 
mentioned  nanowire  devices. 


Figure  1.  7500x  SEM  Micrograph  of  GaN  nanowire  on  Si02 

Substrate 


Traditional  measurement  of  transport  parameters  in 
materials  is  done  electrically^  as  in  diffusion  length 
measurements  performed  through  electron  beam  induced 
current  (EBIC)  experiments.  Electrical  measurements  are 
often  affected  by  contact  related  phenomena  that  require 
potentially  expensive  and  damaging  sample  processing.  In 
view  of  this^  a  contact  free  technique  called  *  transport 
imaging*  was  developed  in  our  laboratory  in  2004  [6] . 

Transport  imaging  employs  an  optical  microscope  (OM)  , 
internal  to  a  scanning  electron  microscope  (SEM) ,  in  order 
to  directly  image  the  spatial  recombination  of  charge 
generated  by  electron  irradiation  at  a  point.  Planar 
contacts  can  be  used  to  create  local  fields  so  that  both 
diffusive  and  drift  behavior  can  be  directly  observed. 

This  has  been  successfully  performed  in  two  dimensional 
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heterostructures  [6,7]  and  in  bulk  materials.  One  goal  of 
this  thesis  is  to  advance  the  technique  by  investigating 
its  initial  applications  to  ID  structures. 

The  SEM-based  approach  also  provides  the  benefit  of 
high  resolution  for  excitation  and  an  ability  to  easily 
excite  wide  bandgap  materials.  The  transport  imaging  method 
makes  it  possible  to  extract  key  materials  parameters,  such 
as  diffusion  lengths  and  minority  carrier  mobility,  from  a 
single  charge  coupled  device  (CCD)  picture. 

B .  MILITARY  RELEVANCE 

Solar  cells  and  integrated  circuits  (ICs)  on 
intelligence  and  communication  satellites  and  laser  target 
designators  and  night  vision  goggles  (NVGs)  in  the 
battlefield  are  just  a  few  of  the  roles  optoelectronic  and 
semiconductor  devices  play  in  the  modern  military.  Their 
importance  cannot  be  overstated,  and  neither  can  the 
importance  of  the  development  of  physical  understanding  and 
processes  to  advance  these  technologies. 

One  of  the  next  steps  in  many  of  these  technologies  is 
the  fabrication  of  nanowire  based  devices.  A  single 
nanowire,  for  example,  can  function  as  a  stand  alone  cavity 
and  gain  medium,  acting  as  a  very  narrow  linewidth  emission 
laser  which  could  have  many  applications,  including  optical 
computing  and  microanalysis  [4]  .  Reduced  lasing  threshold 
currents  (see  Figure  2)  and  diminished  temperature 
dependence  of  the  threshold  current  observed  in  nanowire 
lasers  mean  more  efficient  and  compact  uncooled  laser 
devices  [5]  .  Higher  modulation  factors  in  electro-optical 
modulators  and  new  nanowire  FETs  may  allow  for  faster,  more 
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compact,  efficient  and  inexpensive  optical  integrated 
circuits  (OICs)  and  ICs  [2]  .  As  a  result,  DARPA  and  other 
agencies  that  support  research  and  development  for  military 
optoelectronics  are  investing  significant  resources  in 
nanotechnology . 


Years 

Figure  2 .  Evolution  of  Threshold  Current  for  Semiconductor 

Lasers  (adapted  from  [2]) 

Figure  2  shows  the  impact  of  the  double 
heterostructure  on  threshold  currents  for  semiconductor 
lasers.  The  result  of  the  movement  from  3D  bulk  materials 
to  2D  double  heterostructures  allowed  for  reduction  in 
threshold  currents,  lower  power  consumption  and  a  reduction 
in  physical  dimensions  that  has  put  double  heterostructure 
semiconductor  lasers  in  every  compact  disk  reader  and  many 
other  optoelectronic  devices.  The  military  applications  of 
nanowire  technologies  in  semiconductor  and  optoelectronic 
devices  could  result  in  the  same  widespread  revolution  in 
technology . 
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C .  THESIS  OVERVIEW 

In  this  work  the  transport  imaging  technique  will  be 
expanded  to  examine  spatial  luminescence  from  ID  structures 
in  order  to  develop  a  contact  free  method  of  measuring 
carrier  transport  parameters.  Transport  is  a  key 
application  of  quantum  wires.  The  possibility  of  direct 
measurement  allows  for  more  accurate  transport  parameter 
measurements  and  visualization  of  spatial  variations  along 
the  length  of  the  nanostructure.  Effective  application  of 
the  transport  imaging  method  should  allow  the  determination 
of  values  for  minority  carrier  diffusion  length  and 
potentially,  minority  carrier  mobility.  These  are  important 
for  devices,  such  as  lasers,  that  depend  on  minority 
carrier  injection.  Transport  characterization  in  quantum 
wires  is  still  in  a  relatively  early  stage  and  the 

transport  imaging  technique  would  provide  a  unique  way  to 
extract  important  information. 

Chapter  II  discusses  the  theoretical  and  experimental 
determination  of  the  size  of  the  SEM  beam  as  well  as  the 

dynamics  of  electron  beam-material  interaction  and  carrier 
generation,  critical  to  the  SEM-based  study  of  nanoscale 
structures.  The  one  dimensional  equation  for  steady  state 
minority  carrier  recombination  distribution  solved  for  the 

Gaussian  source  is  derived  and  results  from  numerical 

simulations  based  upon  the  previously  developed  equation 
are  presented. 

In  Chapter  III  the  experimental  apparatus  as  well  as 
the  fabrication  of  multiple  batches  of  wire  structures  on 
Beryllium  (Be)  modulation-doped  p-type  AlGaAs/GaAs/AlGaAs 
heterostructure  using  focused  ion  beam  (FIB)  lithography 
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are  described.  GaN  and  ZnO  quantum  wires  grown  by  chemical 
vapor  deposition  (CVD)  are  also  described  in  detail. 

Chapter  IV  present  the  results  of  the  luminous 
intensity  and  drift  experiments  carried  out  on  the  four 
separate  batches  of  top-down  wire  structures  and  the 
cathodoluminescence  (CL)  spectra  and  results  of  different 
imaging  techniques  on  the  bottom-up  nanowires.  The  final 
chapter  provides  an  analysis  of  the  previously  obtained 
results  and  presents  conclusions  for  both  the  top-down  and 
bottom-up  wire  structures  as  well  as  ideas  for  future 
research . 
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II.  TRANSPORT  IMAGING  THEORY  AND  MODELLING 


A.  INTRODUCTION 

Transport  imaging  utilizes  the  combination  of  a  SEM 
for  high  resolution  charge  generation  and  a  Silicon  CCD  in 
conjunction  with  an  CM  for  the  recording  of  spatially 
resolved  luminescence.  Upon  exposure  to  the  electron  beam, 
a  steady  state  spatial  distribution  of  luminescence 
associated  with  charge  recombination  is  created  which  can 
then  be  recorded  by  the  CCD.  The  CCD  images  can  be  analyzed 
in  order  to  provide  quantitative  measurements  of  local 
minority  carrier  diffusion  and  drift  lengths. 

In  order  to  extract  minority  carrier  transport 
parameters  from  CCD  images,  the  electron  beam-  material 
interaction  and  the  subsequent  generation,  drift,  diffusion 
and  radiative  recombination  of  carriers  in  the  material 
must  be  modeled.  The  SEM  beam's  effects  on  the  sample  can 
be  modeled  as  an  extended  generation  region  within  the 
optically  active  GaAs  layer  of  the  sample.  By  employing 
certain  simplifying  assumptions,  a  model  can  be  developed 
from  the  continuity  equation  which  directly  relates  the 
spatial  luminescence  associated  with  minority  carrier 
recombination  to  the  minority  carrier's  diffusion  and 
drift . 

B.  MODEL  PARAMETERS  AND  ASSUMPTIONS 

1.  Gaussian  Generation  Region  and  Beam  Size 

It  has  been  established  experimentally  that  the  SEM 
beam' s  intensity  follows  a  Gaussian  distribution  to  first 
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order.  Theoretical  calculations  of  the  beam  size  assume 
that  the  beam  diameter,  defined  as  the  Full  Width  at  Half 
Maximum  (FWHM)  of  the  spot  distribution,  is  the  quadrature 
sum  of  the  Gaussian  probe  electron  beam  and  the  various 
aberrations  and  has  a  beam  diameter  of  »^50  nm  at  30  keV  and 
1x10”^°  Amperes  (A)  [8]  . 

In  order  to  model  the  generation  region,  the  goal  was 
to  experimentally  determine  the  beam  diameter  of  our  SEM, 
because  of  potential  deviation  from  theoretical  or 
published  values  during  operating  conditions.  A  Geller 
Microanalytical  Laboratories  Magnification  Reference 
Standard  (MRS-3) ,  pictured  in  Figure  3  below,  was  used  to 
provide  a  high  contrast,  high  resolution  pattern. 


Figure  3.  20x  SEM  Micrograph  of  the  MRS-3  Sample 

Using  the  MRS-3,  the  rate  of  change  of  intensity  along 
sharp  edges  of  the  control  sample  was  correlated  to  the 
diameter  of  the  SEM  beam  experimentally.  The  SEM  beam 
profile  over  a  discontinuity  was  measured  under  a  constant 
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accelerating  voltage  of  30  keV  at  different  probe  currents^ 
as  shown  in  Figure  4 . 


Figure  4.  Probe  Current  Dependence  of  30  keV  SEM  Beam  Profile 

Over  Sample  Discontinuity 

The  beam  intensity  profiles  obtained  experimentally 
were  analyzed  and  an  empirical  value  for  the  SEM' s  beam 
width  was  calculated  as  a  function  of  probe  current.  Values 
for  2o^  where  cr  is  the  standard  deviation  of  the  Gaussian 
distribution^  were  obtained  from  the  positions  of  the  16% 
and  84%  thresholds  for  beam  intensity  and  converted  into 

beamwidths  following  =  2->y2ln^cr  ^  2.35cr  [9]  . 


Probe  Current 

16%-84%  (2a) 

FWHM 

(A) 

(nm) 

(nm) 

1x10'^^ 

34 

40 

3x10'^^ 

39 

46 

1x10'^° 

48 

56 

3x10'^° 

107 

126 

1x10'® 

139 

163 

Table  1.  30  keV  SEM  Beamwidth  Dependence  on  Probe  Current 


9 


The  expected  behavior  as  a  function  of  probe  current 
was  observed  as  increasing  the  probe  current  increased  the 
beam  diameter  and  decreased  the  spatial  resolution  of  the 
SEM.  The  SEM  beam  was  measured  to  be  a  well  focused  beam 
(FWHM  'sfO  nm)  at  low  probe  currents.  A  similar  range  of 
experiments  could  be  carried  out  varying  the  accelerating 
voltage  for  a  given  probe  current. 

Modeling  of  the  beam  size's  relation  to  the  shape  of 
the  luminescent  distribution  shows  a  relative  degree  of 
insensitivity,  especially  when  the  beam  FWHM  <<  Ld,  the 
diffusion  length.  This  is  even  the  case  in  the  ID  diffusion 
scheme,  which  characteristically  shows  increased 
sensitivity  to  generation  region  variations.  Figure  5  shows 
a  comparison  between  ID  generation  regions  of  28  nm  (n=10) 
and  166  nm  (n=60) ,  for  an  Ld  of  3.6  pm. 


Figure  5.  ID  Beamwidth  Comparison,  Ld=3 . 6  pm 
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The  Gaussian  distribution  function  f{x)  =  —j=e”^  ,  where 

\ln: 

2  1 

n  ^  can  then  be  used  to  model  the  generation  region  of 

the  electron  beam  from  the  SEM,  where  the  value  for  n»=30 
corresponds  to  a  30  keV,  1x10"^°  A  electron  beam  of  width  56 
nm.  To  first  order,  in  the  heterostructure  material  of 
interest,  the  generation  region  created  by  the  incident  SEM 
beam  is  approximated  by  the  incident  beam  diameter  on  the 
surface,  since  the  width  of  the  layer  is  only  '»1000  A. 


2.  The  Low  Injection  Limit 

In  order  for  the  spatial  luminescence  information 
recorded  by  the  CCD  to  be  a  direct  indication  of  the 
minority  carrier  distribution  without  resorting  to  more 
complex  models,  the  minority  carriers  must  be  recombining 
into  a  much  larger  and  effectively  constant  number  of 
majority  carriers.  This  is  known  as  the  low  injection 
limit . 

The  number  of  holes  and  electrons  per  unit  volume,  p 
and  n,  can  be  calculated  in  the  optically  active  GaAs 
region  of  the  heterostructure  to  be  ^lO^^cm”^  and  '^lO^cm”^ 
respectively,  at  room  temperature.  Ap  and  An  would  then  be 
the  number  of  holes  and  electrons  per  unit  volume  generated 
by  the  electron  beam  incident  on  the  material.  In  order  for 
the  low  injection  limit  assumption  to  hold  for  the 
modeling,  Ap  p  while  A«  »  n  . 

The  maximum  electron-hole  pair  generation  rate, 
neglecting  electron  beam  energy  loss  due  to  backscattering, 
can  be  calculated  by  G  =  E^IE.,  where  Ei^  is  the  incident 
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beam  energy  and  E.  is  the  ionization  energy  of  an  electron- 
hole  pair  in  GaAs .  For  electron  beam  energies  in  the  5-35 
keV  working  range  of  the  SEM^  E.  is  independent  of  incident 
electron  beam  energy  and  has  been  empirically  determined  to 
be  2.8£’^  [10]^  where  E^  is  the  bandgap  of  GaAs^  -1.42  eV  at 

300  K.  The  total  minority  population  is  then  —Gr  where  I  is 

e 

the  probe  current  and  t  is  the  minority  carrier  lifetime  in 
the  GaAs  heterostructure^  -4  ns. 

The  minority  carrier  concentration  in  the  generation 
region  can  then  be  calculated  by  using  a  hemispherical 
generation  volume  approximation  for  a  30  keV^  10  nA 
electron  beam^  following  the  development  of  Kanaya  and 
Okayama  [11].  This  yields  a  minority  carrier  concentration 
of  -lO^^cm”^.  This  value  is  an  overestimate  of  the  minority 
carrier  concentration  due  to  surface  recombination  and  the 
dimension  of  the  sample  (0.175  pm)  relative  to  that  of  the 
beam  penetration  (2-5  pm)  and  is  still  over  two  orders  of 
magnitude  lower  than  the  effective  majority  carrier 
population  due  to  doping  of  the  GaAs  in  the  samples  of 
interest.  All  experimental  work  will  be  done  below  the  low 
injection  limit  of  a  30  keV^  1  nA  electron  beam. 

Within  the  low  injection  limits  the  model  is  valid^ 
but  some  effects  are  not  physically  observable.  In  the  case 
of  very  short  diffusion  lengths  in  materials  such  as  GaN 
(L-0.1  pm)  that  approach  the  optical  imaging  system's 
resolution  of  -0.4  pm^  the  simulated  effects  cannot  be 
observed.  In  very  large  diffusion  length  materials  such  as 
GalnAs  solar  ceil  heterostructures  (L-130  pm)  it  is  very 
difficult  to  stay  within  the  low  injection  limits  as  the 
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carrier  lifetimes  can  increase  over  three  orders  of 
magnitude,  from  ns  to  ms.  In  addition,  recombination  in 
these  materials  occurs  over  such  a  large  area  and  can  be 
below  the  required  intensity  per  unit  area  to  register 
above  noise  levels  on  the  CCD. 

C.  ID  STEADY  STATE  MODEL 

1 .  The  Non-homogeneous  Second  Order  Differential 
Equation 

The  top-down  wire  structures  described  in  Chapter  III 
vary  from  ID  and  y,z<Lj)  quantum  wires  to  2D  (x,y:s>Lj 

and  z<Lj)  heterostructures  in  their  physical  dimensions. 
The  development  of  steady  state  carrier  recombination 
models  for  all  the  wire  structures  is  required;  however, 
the  model  for  2D  structures  has  been  developed  previously 
[6,7]  .  A  ID  steady  state  equation  is  now  developed, 
starting  with  the  continuity  equation  for  electrons  in  a  p- 
type  material: 

^  =  G.-t/„+l(V.J,)  (1) 

at  q 

where  the  time  rate  of  change  of  electrons  per  unit  volume, 
dn 

— ,  is  expressed  in  terms  of  and  ,  the  generation  and 
dt 

recombination  rates  per  unit  volume,  and  the  divergence  of 
,  the  current  density  vector.  In  the  low  injection  limit 
T,  An 

u^= —  where  An  and  r„  are  the  number  of  mrnorrty  carriers 
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per  unit  volume  and  their  respective  lifetime.  The  current 
density  for  minority  carriers  can  be  defined  as: 


Jn  =  qM„nE  +  qD^Wn 


(2) 


where  ju^  is  the  mobility  of  electrons  in  GaAs^  E  is  the 
applied  electric  field  within  the  material  and  is  the 

diffusion  coefficient  for  electrons.  Substituting  Equation 
(1)  into  Equation  (2)  we  arrive  at: 


^  -  [v  •  {q^„nE  +  qD„Wn)] 

at  r„  q 


(3) 


To  simplify  this  expression,  can  be  related  to  the 


diffusion  Length,  L,  by  D^=  — ,  An^n  since  n  ^1  and,  by 

definition,  for  the  steady  state  there  is  no  time 

dn  ^ 

dependence,  so  —  =  0.  We  also  realize  that  any  gradient  or 

dt 

external  electric  field  will  only  exist  in  one  dimension, 

^  d^  fi 

arbitrarily  x,  such  that  E  =  qEx  and  Equation  (3) 

dx 

can  now  be  written  as: 


0  = 


^dn 

dx 


l}  d^n 
r„  dx^ 


(4) 


Equation  (4)  can  be  further  simplified  by  the  substitution 
S  =  jLi^T^E,  where  S  is  the  drift  length,  and  rearranged  into  a 
non-homogeneous  second  order  differential  equation: 
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(5) 


d^n  S 
dx^  1} 


^  dri^ 
ydxj 


1  G  T 
n  = - 


2 .  General  Solutions  for  a  Point  Source 

The  derivation  of  the  one  dimensional  equation  for 
steady  state  minority  carrier  recombination  distribution  is 

first  solved  for  the  point  source,  such  that  G^=^^S{x) 

where  S{x)  is  the  delta  functional.  Equation  (5)  then 
becomes : 


L^^-S 

dx^ 


^  dn'^ 

ydxj 


■n  =  -gLS{x) 


(6) 


The  Fourier  transform  can  be  taken  of  Equation  (6) 
yielding : 


-L^k^N{k)  -  ikSNik)  -  N{k)  =  -gL 


(7) 


Where  N(k)  is  the  Fourier  transform  of  n.  Equation  (7)  can 
be  re-arranged  to  form: 


W{/c)  =  | 


r 


1 


k  +  CJi  +  C 


(8) 


2  J 


where  the 

constants  c, 

iS 

“r 

and  r  -U  . 

The  operation 

of 

completing 

the  square 

can 

then  be 

performed  on 

the 

denominator 

of  the  expression 

in  Equation 

(8) 
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k^+C,k  +  C^  =  k  +  ^  +C, 

V  2  / 


where  C3  =  • 


arrive  at: 


iS^  +AU 


After  a  change  of  variables  k'=k-\-—^  we 


a  1 

N(k)  =  ^  - ^ . 

zU^’)  +C3 


Equation  (10)  can  then  be  multiplied  by  —  and  inverse 


Fourier  transformed : 


n(x)  = 


{kf  +  C, 


Equation  (11)  can  be  recognized  as  a  Fourier  transform  of 
an  exponential  in  conjunction  with  a  multiplicative 

q. 

- -IX 

exponential  term,  e  ^  ,  yielding  the  one  dimensional 

equation  for  steady  state  minority  carrier  recombination 
distribution  solved  for  the  point  source: 


n{x)  = 


Sx  4s^+AL^ 

re^e 
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3. 


Solution  for  a  Gaussian  Generation  Region 


For  the  case  of  the  source  function,  =  the 

T 

solution  to  the  point  source  can  be  generalized,  as  long  as 

the  function  is  a  general  member  of  a  delta  sequence,  such 

00 

that  lim  [  (x)/(x)c/x  = /(O)  .  The  Gaussian  distribution  function 

-00 

f(x)  =  -^e~”'‘  ,  where  ,  satisfies  these  conditions. 

4^  2c7^ 

Following  the  development  preceding  Equation  (10)  a 
parallel  equation  can  be  derived  which  takes  the  modified 
source  function  into  account: 


mk) 


(13) 


where  F(k)  is  the  Fourier  transform  of  f  (x)  .  The  inverse 
Fourier  transform  of  N(k),  w(x),  can  be  expressed  as  a 

convolution,  that  is  = 


w(x)  - 


m, 

r 


(k')  +c, 


(14) 


Applying  the  convolution  theorem^ 


00  /  00 


Jx' 


(15) 
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and  taking  advantage  of  the  previously  derived  result  in 
Equation  (12)  for  the  result  of  the  Fourier  transform  of 
the  second  term.  Equation  (16)  is  obtained: 


w(x)  =  J 


1  n 


2U 


^/F+4? 


S{x-x') 


+4r 


2L^ 


dx^ 


(16) 


After  rearranging^  Equation  (17)  is  final  form  of  the  one 
dimensional  equation  for  steady  state  minority  carrier 
recombination  distribution  solved  for  the  Gaussian  source: 


w(x) 


uu 

J 


e  e 


S{x-x') 
2L^  , 


+4i^ 


2h 


dx' 


(17) 


D .  MODELING  RESULTS 

1 .  Introduction 

The  steady  state  carrier  recombination  models 
developed  in  the  previous  section  can  be  used  to  calculate 
expected  spatial  intensity  distributions.  MATHCAD  was  used 
to  calculate  the  single  and  double  integrals  numerically, 
employing  a  Romberg  trapezoidal  approximation  with 
constraint  and  convergence  tolerances  of  0.001.  Numerically 
the  integration  over  an  infinite  generation  region  is  not 
possible,  so  a  finite  integration  over  three  standard 
deviations  of  the  Gaussian  generation  region  accounting  for 
99.7%  of  the  total  generation  is  used  as  an  approximation. 
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2. 


Non-Normalized  Diffusion  Modeling 


Non-normalized  diffusion  modeling  of  top-down  wire 
structures  shows  a  substantial  expected  difference  in 
recombination  intensity  between  the  ID  quantum  wires  and  2D 
heterostructures,  assuming  relevant  material  parameters 
stay  constant.  At  its  peak,  the  expected  intensity 
difference  is  approximately  30  times  greater  for  a  ID  drift 
scheme  than  a  2D  drift  scheme,  as  illustrated  in  Figure  6. 
On  the  logarithmic  scale  the  slope  of  the  ID  line  is 
linear,  while  the  function  of  the  2D  line  is  a  zeroth  order 
modified  Bessel  function  of  the  second  kind  whose  slope  is 
approximately  linear  far  from  the  origin. 


Position  (Microns) 

Figure  6.  ID  and  2D  Intensity  Comparison,  Linear  and  Log 

Scales,  Lci=3 . 6  pm,  n=30 


3.  Normalized  Drift  Modeling 

While  the  non-normalized  results  demonstrate  the  clear 
difference  between  ID  and  2D  intensities^  normalization 
allows  for  a  more  effective  comparison  of  the  predicted 
shapes  of  spatial  intensity  distributions.  This  method  of 
comparison  is  useful  because  of  the  absolute  intensity 
fluctuations  due  to  SEM  optics  and  natural  filament 

thermionic  emission  variations  between  experiments.  Figures 
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7  and  8  show  a  change  of  three  orders  of  magnitude  in 
applied  field  on  ID  and  2D  structures  and  the  resulting 
drift  on  both  linear  and  logarithmic  intensity  scales: 


cn 


E 

o 


Figure  7.  ID  Drift  Intensity  Comparison,  Linear  and  Log 

Scales,  Ld=3 . 6  pm,  n=30 


Figure  8.  2D  Drift  Intensity  Comparison,  Linear  and  Log 

Scales,  Ld=3 . 6  pm,  n=30 


4.  Comparison  of  ID  and  2D  Normalized  Drift  Modeling 

In  the  modeling  process  an  arbitrary  amount  of 
electron-hole  pairs  are  created  in  the  generation  region 
and  their  movement  and  recombination  are  simulated.  By 
using  the  same  generation  regions  and  the  same  material 
parameters  in  the  ID  and  2D  models,  only  the  dimensional 
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parameters  affect  the  intensity  and  distribution  of  the 
luminescence . 

In  the  2D  case,  the  generated  charge  drifts  and 
diffuses  into  a  region  (2nrdr)  which  increases  in  direct 
proportional  to  the  distance  traveled.  In  the  ID  case 
however,  there  is  no  change  in  the  size  of  the  region  (dr) 
charge  diffuses  into  with  distance  traveled.  In  practice, 
when  the  carrier  diffusion  lengths  become  much  greater  than 
the  width  of  a  sample  heterostructure,  an  effectively  ID 
drift  and  diffusion  scheme  is  created. 

A  comparison  between  ID  and  2D  normalized  drift 
characteristics  shows  differences  in  the  shape  of  luminous 
intensity  that  increase  with  applied  field.  With  increasing 
drift  ( )  lengths,  proportional  to  increasing  applied 
external  fields,  the  effective  confinement  experienced  by 
the  minority  carries  increases.  Due  to  this  effect, 
transport  in  ID  is  more  sensitive  to  variations  in  material 
parameters.  Figure  9  shows  variation  between  ID  and  2D  low 
and  high  applied  field  drift  conditions  on  both  linear  and 
logarithmic  scales : 


Figure  9.  ID  and  2D  Drift  Intensity  Comparisons,  Linear  and 

Log  Scales,  L;j=  3 .  Gpm,  n=30 
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III.  EXPERIMENTAL  APPARATUS  AND  SAMPLES 


A.  EXPERIMENTAL  APPARATUS 
1 .  The  SEM 

A  JEOL  840A  SEM  with  an  internal  optical  microscope 
(OM)  and  mirror  assembly  allowed  for  SEM,  optical  and  CL 
imaging.  An  OM  objective  connected  to  a  charge  coupled 
device  (CCD)  camera  or  a  mirror  assembly  connected  to  a  CL 
system  was  inserted  by  means  of  retractable  arms  into  the 
SEM  chamber  for  data  collection.  The  SEM  system  uses  a 
constant  flow  liquid  cooled  stage  which  operates  between 
room  temperature  and  about  5K  if  liquid  helium  is  used  as 
the  coolant.  Figure  10  shows  a  schematic  representation  of 
the  SEM  system. 


r 


Monochromator 
and  PMT 


Retractable 
CL  option 


Continu  ous  fl  ow  He 


Figure  10.  Schematic  of  SEM  System 


The  Oxford  Instruments  CL  system  connected  to  the 
retractable  mirror  employs  a  parabolic  mirror,  k  m 

monochromator  and  TE  cooled  GaAs  PMT  as  a  detector.  An 
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Apogee  silicon  CCD  camera  with  a  2184x1472  pixel  array 
serves  as  the  optical  detector  connected  to  the  retractable 
CM. 

The  Apogee  CCD  camera  collects  unfiltered  light  from 
about  300-1100  nm  on  a  15x10  mm  CCD  array.  The  CCD  camera 
is  thermoelectrically  cooled  to  temperatures  below  -10 °C  to 
reduce  noise  during  operation.  Each  pixel  on  the  array  has 
dimensions  of  6.8  pm  x  6.8  pm  and  the  OM  has  a 
magnification  of  approximately  20x.  This  yields  a  CCD  on 
sample  resolution  of  about  0.4  pm/pixel.  The  resolution  for 
incoherent  light  can  be  estimated  as 


0.6U 


NA 


(18) 


where  Ay  is  the  wavelength  dependent  spatial  resolution,  A 
is  the  wavelength  and  NA  is  the  numerical  aperture.  For  the 
room  temperature  peak  emission  of  GaAs  (870  nm)  and  a 
numerical  aperture  of  0.95  (as  a  maximum  value).  Ay  is 
found  to  be  0.56  pm.  Thus  the  CCD  resolution  is  close  to 
the  diffraction  limit  of  870  nm  light  characteristic  of 
GaAs. 


2.  SEM  Modes 

The  SEM  can  operate  in  picture,  line  and  spot  modes. 
Line  and  picture  modes  allow  the  beam  to  raster  in  a 
straight  line  or  over  a  rectangular  area,  while  spot  mode 
focuses  the  beam  down,  as  much  as  possible,  to  a  point.  SEM 
operating  modes  and  magnifications  can  be  chosen  depending 
on  the  specific  generation  region  shape  and  size  desired. 
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In  spot  mode^  the  electron  beam  is  focused  down  to  a 
single  point  whose  actual  size  depends  primarily  on  probe 
current  and  accelerating  voltage  as  discussed  in  Chapter 
II.  Spot  mode  also  displays  a  small  secondary  spot  due  to 
internal  reflection  in  the  OM  optics^  as  seen  in  Figure  11. 


Figure  11.  30  keV,  1x10  A,  1  Second  Spot  on  GaAs 

Heterostructure  (154  pm  x  125  pm) 

Imaging  with  the  SEM  as  well  as  sample  measurements 
can  be  done  in  line  and  picture  modes.  Picture  mode  also 
serves  as  an  ideal  tool  for  fine  focus  adjustments  in 
height  for  the  optical  system  and  can  serve  as  a  method  for 
identifying  important  sample  features  relative  to  a  spot 
when  used  in  conjunction  with  spot  mode  operation  during 
the  same  CCD  exposure. 
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3. 


Data  Extraction 


Diffusion  and  drift  spot  data  in  the  ID  limit  takes  on 
an  exponential  form,  as  described  in  Chapter  II,  which  can 
be  directly  extracted  from  the  value  of  the  slope  in  the 
CCD  pictures  intensity  profiles.  Data  from  the  line  scan 
feature  included  in  the  Micro  CCD  software  package  or  data 
from  an  imported  image  matrix  in  MATLAB  can  be  used  to 
obtain  intensity  profiles  from  CCD  images.  To  reduce  noise 
and  locate  and  extract  luminescence  maxima  from  the  TIF 
files  in  which  CCD  data  is  recorded  a  simple  MATLAB  program 
was  used  (see  Appendix  A) . 

B.  SAMPLE  PROPERTIES  AND  FABRICATION 

1 .  Sample  Heterostructure 

All  top-down  samples  were  fabricated  from  a 
AlGaAs/GaAs/AlGaAs  double  heterostructure  grown  by  solid 
source  molecular  beam  epitaxy  (MBE)  on  a  semi-insulating 
substrate.  Be  was  used  to  modulation  dope  the  AlGaAs 
regions,  making  the  optically  active  GaAs  layer  p-type .  The 
minority  carriers  are  electrons  and  within  the  GaAs  layer 
there  is  an  effective  hole  concentration  of  »^lxlO^^  cm~^ . 
Figure  12  illustrates  the  heterostructure  dimensions. 

Luminescence  from  the  semi-insulating  substrate  is 
negligible  compared  to  that  of  the  epitaxial  GaAs  in  CL 
measurements  due  its  significantly  lower  level  of  optical 
activity  relative  to  the  GaAs . 
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Figure  12.  The  AlGaAs/GaAs/AlGaAs  Heterostructure 

2 .  Focused  Ion  Beam  (FIB)  Lithography  in  Wire 
Structure  Fabrication 

In  order  to  create  the  'top-down'  wire  structures  for 
carrier  diffusion  and  recombination  studies  in  this  thesis, 
one  of  the  methods  used  was  FIB  lithography.  FIB  is  a 
technique  in  the  semiconductor  industry  used  in 
applications  such  as  mask  repair  and  site  specific 
Transmission  Electron  Microscope  (TEM)  sample  preparation. 
Exposing  the  heterostructure  to  the  FIB  served  as  a  means 
of  physically  and  electrically  isolating  the  wire 
structure . 

Inside  the  FIB  system,  a  Ga  liquid  metal  ion  source 
(LMIS)  is  used  to  create  ionized  Ga  ions  which  are  then 
accelerated  through  a  potential  difference  of  10-40  keV  and 
focused  with  electrostatic  lenses.  When  the  high  energy  Ga 
ions  strike  the  sample,  the  surface  is  amorphized  and  atoms 
on  the  surface  are  sputtered,  which  allows  the  sample 
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exposure  to  a  FIB  to  be  analogous  to  milling  the  exposed 
areas . 

The  FIB  system  also  contains  an  SEM  which  operates  at 
an  offset  of  52°  from  the  FIB.  The  sample  is  connected  to  a 
stage  which  can  be  rotated  to  align  it  with  the  electron 
beam  or  the  ion  beam.  Short  of  using  the  ion  beam  to  image, 
however,  the  FIB  etching  must  be  done  while  viewed  from  a 
52°  offset  which  introduces  challenges  to  milling  exact 
structures . 

3.  Chemical  Vapor  Deposition  (CVD)  in  Nanowire 
Fabrication 

The  'bottom-up'  nanowire  structures  studied  in  this 
thesis  were  deposited  by  metal  initiated  metal-organic  CVD 
on  Au  and  Si02  substrates.  The  process  was  carried  out  with 
different  reagents  in  order  to  deposit  both  GaN  and  ZnO 

nanowires  on  the  substrate  at  temperature  of  800~1000  °C 

[14]  .  CVD  is  widely  used  in  the  semiconductor  industry  for 
thin  film  deposition. 

In  CVD  the  substrate  is  placed  on  a  temperature 
controlled  susceptor,  made  from  a  material  resistant  to  the 
metalorganic  compounds.  Nitrogen  gas  is  then  bubbled 
through  a  metalorganic  liquid  and  used  as  a  carrier  gas.  In 
the  case  of  GaN  nanowires,  nitrogen  was  bubbled  through 
Trimethylgallium  (TMG)  to  create  a  metalorganic  vapor  and 
introduced  into  the  sample  chamber  with  ammonia  gas . 

Thermally  evaporated  metals  2-10  nm  thick  on  the 

substrate  such  as  Ni,  Fe  and  Au  are  then  used  as  initiators 

for  vapor-liquid-solid  (VLS)  nanowires.  The  VLS  process 

begins  with  the  dissolution  of  gaseous  metalorganic 
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reactants  into  the  liquid  droplets  of  the  catalyst  metal. 
Nucleation  and  growth  of  singlecrystalline  rods  can  then  be 
achieved  as  the  metalorganic  vapors  combine  in  on  the  seed 
crystal,  allowing  the  growth  of  nanowires  [14] . 

C.  FIRST  SAMPLE  BATCH 

1 .  Batch  Fabrication 

The  fabrication  of  the  first  sample  batch  of  top-down 
structures  began  with  a  methanol  rinse  of  the  fragment  of 
the  AlGaAs/GaAs/AlGaAs  heterostructure  to  be  etched  by  the 
FIB.  In  order  to  create  an  electric  field  within  the  sample 
for  potential  drift  measurements.  Platinum  (Pt)  pads  for 
contacts  were  made  on  the  sample  by  electron  beam 
deposition . 

In  order  to  align  the  wire  structures  with  the  Pt 
pads,  the  sample  was  imaged  by  a  500  pA  ion  beam  for  a 
number  of  seconds  before  it  was  etched  by  the  FIB,  with  8- 
10  minute  exposures  to  a  500  pA  Ga  ion  beam.  The  electron 
beam  deposition  of  Pt  contacts  and  the  subsequent  imaging 
and  etching  by  the  FIB  was  repeated  to  fabricate  all  of  the 
wire  structures. 

The  sample  was  then  rinsed  in  solution  of 
Hydrocholoric  Acid  (HCl)  for  30  seconds  to  remove 
accumulated  Ga  on  its  surfaces  produced  by  the  FIB  etching 
process . 


2 .  Batch  Sample  Dimensions  and  Details 
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The  first  sample  batch  was  fabricated  with  widths 
ranging  from  approximately  1  to  40  pm.  High  resolution  SEM 
micrographs  revealed  variations  in  the  wire  widths  of  up  to 
±100  nm,  as  illustrated  in  Figure  13.  The  irregularities  in 
the  attempted  fabrication  of  uniform  width  structures  also 
brings  first  sample  batch  fabrication  methods  into 
question.  In  Table  2  the  error  column  is  a  combination  of 
both  measurement  uncertainty  and  structure  width  variation. 


Figure  13.  lOOOx  SEM  Micrograph  of  923  nm  Wire  Structure 


Sample 

Wire  Width 

Error 

Error 

Number 

(nm) 

(nm) 

(%) 

1 

39150 

+  100 

0.26 

2 

9977 

+  45 

0.45 

3 

948 

+  35 

3.69 

4 

923 

+  40 

4.24 

Table  2.  First  Batch  Sample  Dimensions 

3 .  First  Batch  Damage  Considerations 

Measurements  on  and  near  the  FIB  etched  wire 
structures  displayed  luminous  intensities  over  1000  times 
lower  than  measurements  taken  a  few  millimeters  away  and 
served  as  a  clear  indication  of  material  damage  associated 
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with  decreased  minority  carrier  lifetime  in  the  vicinity  of 
the  wires.  In  order  to  understand  and  quantify  the  decrease 
in  intensity  in  the  vicinity  of  the  wire  structures,  a 
series  of  CL  spectra  were  taken.  Figure  14  shows  the  four 
FIB  etched  wires,  indicated  by  white  arrows  and  the 
position  of  the  CL  spectra  taken,  indicated  by  the  numbers 
1  to  8 .  The  spectra  data  points  were  evenly  spaced  in 
intervals  of  465  pm. 


Figure  14.  550x  SEM  Micrograph  of  First  Batch  Sample  CL  Points. 
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Figure  15. 


Spatial  Variation  of  CL  Spectra  on  First  Batch 


In  the  vicinity  of  the  wire  structures  CL  luminosity 
decreased  significantly,  and  was  only  measurable  in  the 
areas  where  the  substrate  was  exposed  by  the  FIB  process, 
as  shown  in  Figure  16. 


Figure  16.  SEM  and  Corresponding  CL  Data  for  40  pm  Wire 


Without  the  smaller  bandgap  heavily  doped  GaAs 

material  covering  the  semi-insulating  substrate  in  the  FIB 

etched  regions  it  was  determined  that  observed  luminosity 

came  directly  from  the  substrate  rather  than  the  epitaxial 

GaAs  layer.  Figure  17  shows  broadening  of  the  spectra  and 

shift  in  energy  peaks  typical  of  heavily  doped  materials, 
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when  comparing  the  doped  GaAs  and  the  FIB  etched  well 
substrate  regions  of  the  40  pm  wire  structure. 


Figure  17.  FIB  Etched  Well  Area  and  Epitaxial  Layer  CL  Spectra 

The  data  collected  on  the  first  sample  batch  was  clear 
evidence  of  large  scale  damage  to  the  heterostructure 
material  as  well  as  large  wire  structure  width  variation. 
Exposure  to  the  FIB  imaging  was  enough  to  amorphize  and 
etch  surface  of  the  material  in  the  vicinity  of  the  wire 
structures,  altering  its  optical  properties  significantly, 
more  than  expected  [15] . 

D.  SECOND  SAMPLE  BATCH 

1 .  Batch  Fabrication 

In  order  to  reduce  the  damage  to  the  sample  that 
occurred  in  the  fabrication  of  the  first  batch,  isolating 
the  wire  structure  region  from  FIB  damage  became  a 
priority.  in  order  to  achieve  this,  contacting  on  the 
sample  had  to  be  done  in  a  way  that  did  not  require  FIB 


33 


imaging  of  the  area.  A  200  pm  Silicon  mask  was  affixed 
across  the  selected  AlGaAs/GaAs/AlGaAs  heterostructure 
fragment  with  wax  to  protect  the  area  on  which  the  wires 
would  be  etched  from  electron  beam  deposition  of  Ni/Ge/Au 
contacts  that  followed.  Figure  18  gives  a  schematic 
representation  of  the  area  processing  of  the  sample  batch. 


The  silicon  mask  was  then  physically  removed  and 
methanol  and  xylene  chemical  rinses  were  used  to  clean  the 
area.  Rapid  thermal  annealing  (RTA)  of  the  contacts  was 
then  carried  out  at  430  °C  for  20  seconds  and  the  wire 
structures  made  in  the  non-metalized  region  of  the  sample 
using  8-10  minute  500  pA  FIB  etching.  The  sample  batch  was 
then  sonicated  in  solution  of  HCl,  causing  the  metal 
contacts  to  become  detached  from  the  sample. 

The  contact  areas  were  then  cleaned  with  low  power 
oxygen  plasma  and  the  Si  mask  was  reattached  to  the  sample 
with  wax.  A  second  attempt  at  electron  beam  deposition  of 
Ni/Ge/Au  contacts  also  failed  and  after  sonication  in 
methanol  and  xylene  chemical  rinses  the  sample  was  cut  with 

a  wire  saw,  to  electrically  isolate  each  wire  structure,  as 
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show  in  Figure  18.  Pressed  In  contacts  were  added  to 
individual  samples  as  required  for  application  of  external 
bias . 


2 .  Batch  Sample  Dimensions  and  Details 

The  second  sample  batch  was  fabricated  with  wire 
structure  widths  which  ranged  from  about  0.7  to  24  pm.  No 
visual  evidence  of  amorphization  on  the  sample  surface  due 
to  exposure  to  the  FIB  was  present,  but  partially  deposited 
contacts  which  failed  to  adhere  to  the  sample  surface  are 
clearly  visible  in  Figure  19. 


Sample 

Number 

Wire  Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

23630 

+  85 

0.36 

2 

13890 

+  50 

0.36 

3 

1725 

+  20 

1.16 

4 

746 

+  35 

4 . 69 

Table  3.  Second  Batch  Sample  Dimensions 
Figure  19.  lOOX  SEM  Micrograph  of  Contacted  Second  Sample  Batch 


A  significant  increase  in  the  consistency  of  the  wire 
widths  was  also  observed,  as  shown  in  Figure  20. 
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Figure  2  0  . 


7500x  SEM  Micrographs  of  1.7  pm  and  14  pm  Structures 


3 .  Second  Batch  Damage  Considerations 

In  the  second  sample  batch  optical  measurements  near 
the  FIB  etched  wire  structures  displayed  luminous 
intensities  over  10  times  higher  than  measurements  taken  in 
the  same  regions  of  the  first  sample  batch.  While  less 
damage  was  done  to  the  samples^  the  absolute  intensity  in 
these  regions  was  still  almost  100  times  lower  than  that  of 
the  control  sample. 

The  700  nm  wire  demonstrated  both  light  and  dark 
regions  under  SEM  and  CL  measurements  and  visible  variation 
in  its  width  and  structure.  All  wire  structures  also 
demonstrated  surface  irregularities  in  the  regions 
surrounding  the  FIB  milling.  The  second  batch  samples  were 
exposed  to  multiple  RTA  and  chemical  processes  and  the 
observed  data  indicate  that  these  processes  could  have 
damaged  or  degraded  the  materials.  Figure  21  shows  the 
surface  irregularities  in  the  vicinity  of  the  700  nm  wire 
as  well  as  the  irregularities  within  it. 


Figure  21.  550x  and  12000x  SEM  Micrographs  of  746  nm  Wire 

Damage 

E.  THIRD  SAMPLE  BATCH 

1 .  Batch  Fabrication  Details 

The  third  sample  batch  was  fabricated  with  a  focus  on 
avoiding  both  overprocessing  and  direct  exposure  to  the  FIB 
etching,  the  presumed  causes  of  the  decreased  luminescence 
in  previous  batches .  A  heterostructure  sample  was  cut  into 
five  parts  and  in  order  to  further  test  damage  theories  a 
pressed  In  protective  cover  was  placed  on  one  of  the 
samples.  All  of  the  samples  were  then  put  in  the  FIB 
chamber . 

A  50  pm  X  100  pm  area  was  FIB  etched  at  500  pA  for  0.2 
seconds  and  another  23  pm  x  33  pm  area  was  FIB  etched  at 
500  pA  for  2  minutes  on  the  In  protected  sample.  Two  =^0.5 
pm  wire  structures  were  also  etched  on  this  sample  using 

both  an  8  minute  500  pA  FIB  exposure  and  40  minute  30  pA 

FIB  exposure. 

Wire  structures  were  then  FIB  etched  into  three  of  the 
remaining  four  samples,  with  the  fourth  serving  as  a 

control  sample.  The  samples  were  then  rinsed  in  HCl 

solution  for  30  seconds  and  pressed  In  contacts  were  added 


to  individual  samples  as  required  for  application  of 
external  bias . 

2 .  Batch  Sample  Dimensions 

The  third  sample  batch  was  fabricated  with  one  larger 
pm)  and  three  smaller  (^0.5  pm)  wire  structures.  A 
significant  increase  in  the  consistency  of  the  wire  widths 
was  also  observed  over  the  last  two  batches,  and  one  300  nm 
wire  was  constructed.  Average  sample  dimensions  were 
acquired  through  line  scans  from  the  SEM  and  the  degree  of 
uncertainty  in  the  measurements  correlates  directly  to  the 
intensity  contrast,  however,  in  the  case  of  the  300  nm 
wire,  the  resolution  limitations  of  the  SEM  rather  than  the 
wire  quality  determined  the  errors  in  measurement. 


Sample 

Number 

Wire  Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

6529 

+  35 

0.54 

2 

489 

+  15 

3.07 

3 

571 

+  20 

3.50 

4 

300 

+  30 

10.0 

Table  4.  Third  Batch  Sample  Dimensions 


Figure  22.  lOOx  SEM  Micrograph  of  Third  Batch  Samples 
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F.  FOURTH  SAMPLE  BATCH 

1 .  Batch  Fabrication 

The  fourth  batch  of  wire  structures  was  etched  with  a 
FIB  directly  into  the  already  contacted  and  functional 
sample  that  was  used  as  a  control  in  the  third  batch  of 
experiments.  This  was  done  to  avoid  contacting  difficulties 
present  in  previous  samples  and  to  have  extensive  pre-FIB 
exposure  drift  and  intensity  measurements  for  comparison 
with  the  fourth  batch  samples . 

The  use  of  box  templates  in  the  FIB  patterning 
software  allowed  for  a  more  accurate  FIB  pattern  on  the 
sample,  eliminating  the  width  variations  occasionally 
produced  by  line  template  patterning  used  previously.  The 
sample  was  etched  with  10  minute,  500  pA  series  exposures 
of  the  FIB. 

2 .  Batch  Sample  Dimensions 

A  series  of  four  wires  were  fabricated  on  the  fourth 
sample  batch,  ranging  in  width  from  approximately  0.8  to  9 
pm.  Figure  23  shows  a  high  definition  lOOOx  SEM  micrograph 
of  the  fourth  batch  samples  etched  into  the  GaAs  double 
heterostructure.  The  vertical  dimension  of  the  wires  is 
1 5  3 . 5  pm . 
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Figure  23. 


lOOOx  SEM  Micrograph  of  Fourth  Batch  Sample 


Average  sample  dimensions^  shown  in  Table  5^  were 
acquired  through  line  scans  from  the  SEM  and  the  degree  of 
uncertainty  in  the  measurements  correlates  directly  to 
resolution  limitations  of  the  SEM. 


Sample 

Number 

Wire  Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

9258 

+  20 

0.36 

2 

4850 

+  25 

0.36 

3 

1632 

+  20 

1.16 

4 

792 

+  20 

4 . 69 

Table  5.  Fourth  Batch  Sample  Dimensions 
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IV.  EXPERIMENTAL  RESULTS 


A.  INTRODUCTION 

The  results  of  the  characterization  work  can  be  broken 
down  into  two  main  categories,  intensity  and  drift 
measurements  on  the  top-down  GaAs  wire  structures  and  the 
initial  observation  of  luminescence  from  the  bottom-up  GaN 
and  ZnO  nanowire  structures.  The  top-down  wire  structure 
results  are  further  divided  into  their  respective  batches, 
while  the  bottom-up  structures  are  divided  into  SEM  and  CL 
imaging,  CL  spectra  and  CCD  imaging. 


B.  TOP-DOWN  WIRE  STRUCTURES 

1 .  First  and  Second  Sample  Batches 

The  first  and  second  batches  represented  initial 
attempts  at  production  of  samples  displaying  ID 
confinement.  As  discussed  in  the  previous  chapter, 
significant  damage  was  observed  after  FIB  etching  and 
recorded  via  luminescence  decreases  in  both  the  first  and 
second  batch  samples.  This  was  initially  attributed  to  the 
multi-step  processing  that  the  samples  underwent. 

Spot  intensity  measurements  were  done  on  the  second 
batch  samples  and  an  overall  pattern  of  decreasing  average 
intensity  was  observed  as  the  wire  structure  dimensions 
decreased.  Average  spot  intensities  are  shown  for  the  four 
different  second  batch  sample  wire  structures  under  the 
same  electron  beam  excitation  of  1x10”^°  A  at  30  keV  in 
Figure  24 . 
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Log(Wire  Width)  (Microns) 

Figure  24.  Comparison  of  Average  Second  Batch  Wire  Structure 

Intensity 


2 .  Third  Sample  Batch 


The  third  sample  batch  provided  a  quantitative 
measurement  of  the  FIB  damage  to  the  samples,  a  comparison 
of  wires  etched  at  different  FIB  currents  and  initial 
contacting  and  drift  measurements.  The  third  batch  control 
sample  was  successfully  contacted  and  drift  measurements 
were  carried  out  on  it,  but  those  results  are  presented  in 
the  next  section  due  to  their  relevance  to  batch  four. 


The  third  batch  samples  were  contacted  using  pressed 
In  contacts  and  the  resistances  of  the  contacted  samples 
were  measured  using  a  fine  tipped  probe  station.  Resistance 
values  were  far  below  the  control  sample  scaled  values,  as 
displayed  in  Table  6.  The  significantly  lower  than 
predicted  resistances  and  the  inability  to  observe  drift  in 
the  third  batch  samples  led  to  concern  for  the  quality  of 
electrical  contacts  on  the  samples  and  raised  the 
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possibility  of  parallel  conduction  paths  due  to  surface 
damage . 


Sample 

Control 

2  4  ym 

14  ym 

1 . 4  ym 

0 . 7  ym 

Measured  Resistance 

(kQ) 

0.3 

7 

9.4 

41 . 6 

120 

Predicted  Resistance 

(kQ) 

N/A 

18.5 

31 

251 

592 

Table  6.  Third  Batch  Contacted  Sample  Resistance 


Using  an  In  protective  layer  on  one  of  the  samples 
allowed  for  three  regions  of  different  FIB  exposure  to  be 
tested.  A  50  ym  x  100  ym  area  was  FIB  etched  at  500  pA  for 
0.2  seconds  and  a  second  23  ym  x  33  ym  area  was  FIB  etched 
at  500  pA  for  2  minutes.  The  CCD  image  in  Figure  25  shows 
the  intensity  contrast  between  the  0.2  and  120  second  FIB 
exposed  regions  (light  and  dark  respectively) .  The  average 
intensity  results  on  the  In  protected  and  FIB  exposed 
regions  of  the  sample  are  displayed  in  Table  7. 


Figure  25.  2000x,  30  keV,  1x10  A,  5  Second  SEM  Picture  Mode 

CCD  Exposure 
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Region  of  Interest 

Average  Intensity  (CCD  Counts) 

In  Protected  Region 

7201 

0.2  Sec,  500  pA  Exposure 

2391 

120  Sec,  500  pA  Exposure 

1994 

Table  7.  Average  Maximum  0.25  Second  Spot  Intensities,  30  keV,  1x10  A 

Third  batch  samples  also  served  to  examine  the 
intensity  differences  of  wires  etched  at  different  FIB 
currents.  Table  8  displays  the  average  intensity  from  the 
489  nm  8  minute  500  pA  FIB  exposure  and  the  571  nm  40 
minute  30  pA  FIB  exposure  wire  structures.  This  was 
evidence  that  the  exposure  amperage  in  the  FIB  did  not 
change  the  luminescence  of  the  wire  structures 
significantly . 


Region  of  Interest 

Average  Intensity  (CCD  Counts) 

489  nm^  500  pA  Exposure 

544 

571  nm^  30  pA  Exposure 

584 

Table  8.  Average  Maximum  1  Second  Spot  Intensities,  30  keV,  1x10  A 

3 .  Fourth  Sample  Batch 

In  order  to  control  as  many  factors  as  possible  and 
potentially  observe  carrier  drift  in  the  samples^  the 
fourth  batch  of  wires  was  etched  with  a  FIB  directly  into 
the  already  contacted  and  functional  sample  that  was  used 
as  a  control  in  the  third  batch  of  experiments.  Fourth 
batch  sample  wire  structures  displayed  higher  luminous 
intensities  and  greatly  decreased  variations  in  width  due 
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to  the  refinement  of  the  FIB  etching  process  employed  to 
fabricate  them. 

Extensive  spot  intensity  measurements  were  done  on  the 
fourth  batch  samples.  Consistent  with  the  second  batch 
observations,  as  the  wire  structure  dimensions  decreased, 
an  overall  pattern  of  decreasing  average  intensity  was 
observed.  The  average  intensity  observed  on  the  widest 
structure  is  a  factor  of  less  than  the  initial  (pre-FIB 

processed)  material.  The  average  spot  intensities  in  the 
four  fourth  batch  sample  wire  structures  are  shown  in 
Figure  2  6 . 


400 


Log(Wire  Width)  (Microns) 


Figure  26. 


Comparison  of  Average  Fourth  Batch  Wire  Structure 

Intensity 


In  order  to  analyze  drift  in  the  fourth  batch  wire 

structures^  it  is  important  to  first  examine  drift  data 

from  the  contacted  heterostructure  used  as  the  control 

sample  in  the  third  batch  experiments.  Figure  27  shows  data 
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obtained  from  the  sample  before  it  was  exposed  to  the  FIB 
at  fields  of  zero  and  ^0.01  V/pm  and  Figure  28  shows  line 
scans  of  the  same  unexposed  sample. 


Figure  27.  30  keV,  1x10  A,  0.35  Second  Spot  on  GaAs 

Heterostructure,  0  and  20V  bias  (170  pm  x  160  pm) 


20  40  60  so  100  120  60  80  100  120  140 


Position  (Microns)  Position  (Microns) 

Figure  28.  30  keV,  1x10”^°  A,  0.35  Second  Spot  Linescans  on  GaAs 

Heterostructure,  0  and  ±0.01  V/pm  bias 

Using  the  same  contacts,  potential  differences  of  up 
to  20  volts  were  applied  to  the  sample  after  the  FIB 
etching  of  wire  structures.  A  sample  resistance  of  908  Q 
was  measured  and  there  was  optical  evidence  of  an  applied 
field.  Changes  in  spot  intensity  from  152  counts  per  second 
at  no  bias  to  33  counts  per  second  at  20V  bias  also  serve 
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as  a  clear  indication  that  field  was  applied  to  the  sample. 
Assuming  that  the  voltage  dropped  at  the  contacts  remained 
similar,  an  approximate  potential  difference  of  =«0.13  V/pm 
was  applied  during  experiments.  Figure  29  shows  two  CCD 
images  taken  at  different  bias  on  the  smallest  ('sQ.S  pm) 
wire  structure  and  Figure  30  shows  line  scans  across  the 
same  data.  These  results  are  further  discussed  and  analyzed 
in  Chapter  V. 


Figure  29.  30  keV,  1x10  A,  20  Second  Spot  on  0.8  pm  wire 

structure,  0  and  20V  bias  (136  pm  x  117  pm) 


wire  structure,  0  and  ±0.13  V/pm  bias 
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C.  BOTTOM-UP  WIRE  STRUCTURES 

1 .  SEM  and  CL  Micrographs 

The  first  step  in  expanding  the  transport  imaging 
technique  to  nanowires  involved  locating  the  nanowires  on 
their  respective  substrates  and  recording  SEM  and  CL 
micrographs  to  determine  the  spectra  and  magnitude  of 
luminous  emission.  It  is  important  to  note  that  the  CL 
system  collects  all  light  emitted  from  any  point  in  the 
image  and  maps  it  to  the  point  on  the  image  that  is  being 
exposed  to  the  electron  beam,  so  it  is  used  to  provide  a 
general  idea  of  the  panchromatic  intensity  of  the  sample  or 
the  intensity  of  the  sample  at  a  selected  wavelength,  but 
its  spatial  information  assumes  that  charge  generation  and 
luminescence  arise  from  the  same  point. 

ZnO  and  GaN  nanowires  were  deposited  on  Gold  and  Si02 
substrates  but  could  only  be  located  on  the  Si02  substrate, 
possibly  due  to  chemical  interaction  and  decomposition  of 
the  nanowires  with  the  Au  substrate.  ZnO  and  GaN  nanowires 
and  nanowire  groups  were  observed  to  be  bright  on  SEM 
micrographs  mainly  due  to  charging  because  of  their 
electrical  isolation  on  the  Si02  substrate.  Figures  31  and 
32  show  SEM  micrographs  of  ZnO  and  GaN  nanowire  groups  at 
2000x  and  isolated  nanowires  at  a  magnification  of  4000x  on 
the  Si02  substrate.  The  measured  diameters  of  10-200  nm  and 
aspect  ratios  of  10-100  of  the  nanowires  were  consistent 
with  expected  values  [17,18,19]  for  both  GaN  and  ZnO. 
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Figure  31.  2000x  and  4000x  SEM  Micrographs  of  ZnO  Nanowires,  15 

keV,  1x10'^°  A 


Figure  32.  2000x  and  4000x  SEM  Micrographs  of  GaN  Nanowires,  15 

keV,  1x10"^°  A 

After  the  initial  SEM  micrographs  were  taken,  the  CL 
system  was  activated  in  parallel  with  the  SEM,  allowing  for 
simultaneous  SEM  and  CL  imaging.  The  CL  system  was  set  to 
panchromatic  mode  in  order  to  achieve  maximum  intensity  for 
the  micrographs.  The  results  for  the  ZnO  nanowires  are 
displayed  in  Figure  33. 
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2000x  SEM  (left)  and  CL  (right)  Micrographs  of  ZnO 
Nanowires,  15  keV,  1x10"^°  A 


nanowires  in  panchromatic  mode^  as  shown  in  Figure  34.  The 
panchromatic  CL  image  was  nearly  featureless  in  the  GaN 


case  at  magnifications  of  2000x  due  to  strong  substrate 


emissions  in  the  background 


however  some  individually 


luminescent  wires  were  visible 


2000x  SEM  (left)  and  CL  (right)  Micrographs  of  GaN 
Nanowires,  15  keV,  1x10"^°  A 


luminescent 


regions  on  the  GaN  nanowire  sample  at  higher  magnification 


led  to  successful  CL  imaging  on  that  sample.  Individual  GaN 


nanowires  are 


Figure  35.  lOOOOx  SEM  (left)  and  CL  (right)  Micrographs  of  GaN 

Nanowires,  15  keV,  1x10"^°  A 

2 .  Cathodoliaminescence  Spectra 

CL  spectra  were  taken  on  the  ZnO  and  GaN  nanowire  and 
substrate  regions  in  order  to  determine  the  wavelength 
distribution  of  the  panchromatic  data  taken  in  the  CL 
micrographs.  CL  spectra  were  taken  between  300  nm  and  700 
nm  with  steps  of  1  nm^  to  identify  characteristic  nanowire 
emission  and  distinguish  it  from  the  substrate  emission. 

ZnO  nanowires  displayed  a  very  strong  peak  at  about 
378  nm  (3.29  eV) ^  visible  in  Figure  36. 


Figure  36.  CL  Spectra  of  ZnO  Nanowires  and  Si02  substrate,  15 

keV,  1x10"^°  A 


51 


The  GaN  nanowire  spectra  had  to  be  directly  compared 
to  substrate  spectra  in  order  to  clearly  indentify  an 
emission  peak.  In  Figure  37,  a  small  deviation  from  the 
substrate  spectra  can  be  seen  in  the  GaN  nanowire  spectra 
in  the  vicinity  of  350-450  nm  and  corresponds  to  a  peak  at 
about  380  nm  (3.27  eV)  . 
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Figure  37.  CL  Spectra  of  GaN  Nanowires  and  Si02  substrate,  15 

keV,  1x10"^°  A 

The  observed  spectral  peaks  associated  with  the 
nanowires  are  in  agreement  with  expected  values  of  ZnO  and 
GaN  bandgaps  of  «^3.3  eV  and  «^3.4  eV  respectively  at  300K 
[16,17]  . 


3 .  CCD  Imaging  of  Nanowires 

The  data  gathered  during  CL  clearly  indicate  that  the 
nanowires  luminesce;  however,  the  CL  system  employs  a  high 
voltage  photomultiplier  while  the  CCD  system  records 
directly  emitted  intensity  without  any  optical 
amplification.  In  order  for  the  transport  imaging  technique 
to  be  employed  on  nanowires  they  must  be  visible  on  the  CCD 
system . 
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The  final  results  taken  on  the  nanowire  samples  were  a 
series  of  unfiltered  CCD  images  in  which  GaN  and  ZnO 
nanowires'  spatial  luminescence  was  resolved  with  the  OM 
and  CCD  system  used  for  transport  imaging.  In  Figure  38  a 
group  of  GaN  nanowires  are  imaged  with  the  CCD,  in  Figure 
39  two  single  ZnO  nanowires  are  imaged. 


Figure  38.  2500x,  15  keV,  1x10  A,  120  Second  Image  of  GaN 

Nanowires 
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Figure  39. 


2000x,  30  keV,  1x10  A,  120  Second  Image  of  ZnO 
Nanowires 
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V. 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER 


RESEARCH 

A.  SUMMARY  AND  CONCLUSIONS 

1 .  Top-Down  Wire  Structures 

Throughout  the  research  on  top-down  wire  structures 
four  separate  sample  batches  were  fabricated  and  various 
intensity  and  drift  measurements  were  made  on  the  samples. 
Extensive  experience  with  the  FIB  system  and  samples 
allowed  for  a  progressive  refinement  of  the  FIB  etching 
process  employed  to  fabricate  the  structures,  however, 
there  was  a  significant  decrease  in  average  luminous 
intensity  in  all  samples  that  were  exposed  to  the  FIB.  Many 
wire  structures  displayed  evidence  of  electrical  contacting 
problems  and  even  those  which  were  contacted  effectively 
did  not  display  measurable  drift. 

Initial  FIB  etching  resulted  in  average  decreases  in 
luminous  intensity  of  over  a  factor  of  1000  from  the 
unexposed  AlGaAs/GaAs/AlGaAs  double  heterostructure  and 
wire  structures  of  poorly  defined  edges  and  widths. 
Reduction  in  the  processing  of  the  sample  and  experience 
with  the  FIB  system  and  the  patterning  software  yielded 
significant  increases  in  the  wire  structures'  quality  and 
reduced  the  average  decreases  in  luminous  intensity  to  a 
factor  of  65. 

Although  the  latter  sample  batches'  loss  of  intensity 
and  structural  variation  was  significantly  lower  than 
earlier  batches,  all  the  top-down  wire  structure  batches 
demonstrated  a  consistent  pattern:  Figure  40  shows  that  in 
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all  measured  samples^  as  the  wire  structure  dimensions 
decreased  so  did  the  average  intensity. 
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Figure  4  0 . 


Comparison  of  Average  Second  and  Fourth  Batch  Wire 
Structure  Intensities 


Changes  in  electrical  properties  of  samples  due  to  FIB 
exposure  have  been  recorded  in  the  past  and  characterized 
as  ion-induced  surface  defects  and  passivation  of  dopants 
due  to  surface  amorphization  [20^21],  ion  induced  surface 
defects  could  serve  as  extra  scattering  centers  in  the 
vicinity  of  the  FIB  etching  which  could  change  the  non- 
radiative  lifetime  of  the  carriers  within  the  material  and 
subsequently^  to  first  order^  the  overall  lifetime  and 
luminous  intensity  observed.  This  can  be  modeled  as  an 
effective  decrease  in  diffusion  lengthy  following: 


I. 


(18) 


1 

T 


(19) 


where  p  is  the  mobility,  i  is  the  lifetime,  ir  is  the 
radiative  lifetime  and  inr  is  the  non-radiative  lifetime 
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associated  with  an  increase  in  scattering  centers  and  point 
defects  in  the  material. 

Examination  of  the  fourth  top-down  wire  structure 
batch  allowed  a  direct  comparison  between  previously 
undamaged  material  which  displayed  drift  and  the  same 
material  after  the  FIB  process.  The  unexposed  sample 
displayed,  on  average  '»26,  000  counts  per  second  on  the  CCD, 
while  the  same  regions  after  indirect  exposure  to  the  FIB 
had  intensities  averaging  ^400  counts  per  second,  a  factor 
of  65  decrease  in  intensity. 

Employing  Equations  (18)  and  (19)  a  new  Ld  of  0.44  pm 
can  be  calculated  under  the  assumption  that  to  first  order, 
changes  in  intensity  can  be  interpreted  as  changes  in 
lifetime.  Figure  41  shows  a  drift  model  using  these 
modified  parameters  compared  with  the  experimental  results 
for  smallest  (=^0.8  pm)  fourth  batch  wire  structure. 


Figure  41.  0  and  0.13  V/pm  bias  Ld=0.44  pm,  n=30  Simulation 

(left)  and  30  keV,  1x10”^°  A,  20  Second  Spot  Linescan 

While  there  was  a  clear  broadening  of  the  normalized 
spot  under  applied  bias  in  both  directions,  there  was  no 
reliable  evidence  of  directional  drift  as  illustrated  in 
the  modeling.  The  main  limitation  to  drift  measurements  on 
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this  scale  is  the  lack  of  optical  resolution  on  the  CCD 
system^  as  even  modeled  data  displayed  directional  drift 
>1.5  less  than  4  pixels  on  the  CCD  image. 

The  change  in  effective  drift  length  due  to  a  change 
in  lifetime  of  the  minority  carriers  produces  a  model  which 
predicts  drift  below  the  resolution  of  our  optical  system^ 
which  is  consistent  with  the  observed  data.  This  hypothesis 
also  effectively  predicts  the  observed  trend  of  decreasing 
intensity  with  decreasing  wire  dimensions  because  if  a 
constant  area  around  the  FIB  etched  areas  is  damaged  then 
the  smaller  wire  structures  would  display  a  decreased 
luminous  intensity  as  a  larger  fraction  of  their  area  would 
have  been  damaged. 

Although  the  methods  of  FIB  etching  sample  fabrication 
improved  dramatically  throughout  the  course  of  research  and 
the  loss  of  average  luminous  intensity  in  the  samples  was 
significantly  reduced^  we  can  conclude  that^  on  the  current 
apparatus^  until  further  work  has  been  done  to  characterize 
FIB  damage^  FIB  etched  top-down  fabrication  on  GaAs  double 
heterostructure  material  is  not  a  viable  method  of  imaging 
ID  confinement  and  employing  the  transport  imaging 
technique.  Throughout  the  sample  testing^  reduction  in  the 
average  luminous  intensity  and  the  correlated  reduction  in 
diffusion  length  yielded  changes  below  the  optical 
resolution  thresholds  for  our  experimental  apparatus. 


2 .  Bottom-Up  Wire  Structures 


ZnO 

and 

GaN  nanowires 

on  Si02  substrates 

were 

identified 

and 

imaged 

using  the 

SEM,  panchromatic  CL 

and 

the  CCD 

camera . 

Nanowire 

luminescence  was 

also 
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characterized  by  taking  CL  spectra  of  the  nanowires  and  the 
Si02  substrate.  ZnO  and  GaN  nanowires  on  Au  substrates  were 
not  observable . 

Both  the  GaN  and  the  ZnO  wires  showed  spectral  maxima 
that  occurred  in  the  vicinity  of  380  nm,  at  a  higher  energy 
than  most  of  the  broadband  spectra  of  the  Si02  substrate. 
Short  wavepass  filters  could  easily  be  employed  to  reduce 
the  amount  of  background  light  in  the  system  during 
measurements  and  further  decrease  the  substrate  intensity 
(noise)  in  CCD  pictures  in  future  experiments. 

Further  research  into  the  Au  substrates'  chemical 
interaction  with  the  GaN  and  ZnO  nanowires  to  determine  the 
causes  of  nanowire  disintegration  or  lack  of  luminescence 
could  also  be  pursued,  however  this  is  less  critical 
because  future  wire  structures  will  be  on  Si02  substrates. 
We  can  conclude  that  luminescence  of  the  quantum  wires  can 
be  imaged  and  that  their  emission  spectra  are  consistent 
with  expected  values.  This  data  serves  as  a  strong 
indication  of  the  feasibility  of  expanding  the  transport 
imaging  process  to  ID  bottom-up  nanowires. 

B.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

1 .  Advances  in  Bottom-Up  Nanowire  Structures 

The  fabrication  of  contacted  nanowires  is  the  next 
step  in  the  application  of  the  transport  imaging  method  to 
ID  nanowires.  Currently  nanowires  are  grown  on  the 
substrate  surface  by  CVD  and  remain  uncontacted.  Research 
groups  fabricating  single  nanowire  devices  have  reported 
the  ability  to  move  individual  nanowires  and  contact  them 
[4,18,19]  . 
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Contacted  nanowires  would  allow  for  initial 
qualitative  observation  of  drift  in  ID  structures,  and 
coupled  with  increased  system  optical  resolution  could 
yield  the  first  quantitative  drift  measurements.  High 
resolution  drift  measurements  could  then  be  used  to  extract 
transport  parameters  directly  from  CCD  images  of  the 
nanowires . 


2 .  Increased  System  Resolution  for  Transport  Imaging 

In  this  work,  the  use  of  an  optical  microscope, 
internal  to  an  SEM,  to  directly  image  the  spatial 
recombination  of  charge  generated  by  electron  irradiation 
in  nanowires  has  been  demonstrated.  However,  the  current 
system  resolution  is  '«400  nm/pixel  and  even  in  these 
preliminary  examinations  of  nanowire  structures,  the 
apparatus  has  been  pushed  near  the  limits  of  its  optical 
resolution . 

A  new  National  Science  Foundation  grant  has  been 
obtained  to  advance  the  technique  to  higher  resolution 
using  near  field  optical  techniques.  A  near  field 

microscope,  with  special  open  architecture  for  use  inside  a 
SEM,  could  increase  resolutions  of  '^50-100  nm,  and  has  been 
ordered  from  Nanonics  Inc.  and  should  be  installed  by 
summer  2006. 

Increased  system  resolution  would  be  a  very  important 
step  towards  quantitative  transport  imaging  measurements  in 
nanowire  structures  as  well  as  in  materials  which  display 
very  short  diffusion  lengths  such  as  GaN.  System 
resolutions  of  50-100  nm  would  also  allow  for  empirical 
verification  of  the  reduced  diffusion  length  hypothesis  in 
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FIB  etched  top-down  structures  and  for  imaging  of  ID 
confinement  in  bottom-up  nanowire  structures. 

3 .  Transport  Imaging  as  a  Measurement  of  FIB  Damage 

Most  analyses  of  FIB  damage  use  MonteCarlo  simulations 
models  and  analysis  of  damage  layer  thickness  and  damage 
processes  by  TEM  [20] .  TEM  studies  have  also  shown  that  the 
damage  layers  are  amorphized  by  the  FIB,  but  results 
gathered  from  the  top-down  wire  structures  indicate  that 
FIB  etching  may  produce  electrical  damage  that  is  more 
extensive  than  structural  damage  visible  by  TEM. 

An  in  depth  exploration  of  FIB  damage  by  means  of 
transport  imaging  could  provide  a  very  sensitive 
measurement  of  sample  damage  at  an  electronic  or  point 
defect  scale.  Quantitative  studies  would  have  to  be 
undertaken  in  different  materials  and  different  FIB 
exposure  times  and  energies. 

More  research  and  experimentation  would  also  have  to 
be  done  to  determine  if  the  observed  decreases  in  luminous 
intensity  do  in  fact  correspond  to  decreases  in  lifetime 
within  the  material.  Confirmation  of  the  hypotheses 
presented  in  this  paper  could  involve  independent  Time 
Resolved  Photoluminescence  (TRPL) ,  in  which  the  sample  is 
excited  by  a  picosecond  pulsed  laser  and  radiative  carrier 
lifetime  can  be  accurately  measured. 
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APPENDIX  A:  DATA  EXTRACTION  PROGRAM 


1  %Data.m  Made  by  Stephen  winchell 

2  %Last  Modified,  Feb  3,  2006 

3  %Function  to  extract  the  Row  and  Col  Max  from  CCD  Image 

4  %Maxima  marked  by  cross  for  visual  error-check 

5  %Median  noise  filtering  and  noise  subtraction  from  max  value 

6 

7  function  [colmax,  rowmax] =data (pic) ; 

8 

9  %input  image  file  to  matlab  and  convert  it  to  double  to  work  with 
later 

10  I  =  imread(pic); 

11  I  =  double  ( I )  ; 

12  imagesc  ( I ) 

13 

14  %do  a  median  noise  filtering  on  the  image 

15  J  =  medfilt2 (I) ; 

16  imagesc ( J) 

17 

18  %take  off  the  noise  floor  by  subtracting  the  average  of  the 
smallest  sum  of  a  row  or  col 

19 

20  [rows,  cols]  =  size(J); 

21  colsum  =  sum(J); 

22  rowsum  =  sum(J'); 

23 

24  if  (min ( colsum) /cols )  <=  (min (rowsum) /rows ) 

25  noisefloor  =  min (colsum) /cols ; 

26  else 

27  noisefloor  =  min ( rowsum) /rows ; 

28  end 

29 

30  %find  the  maxima  and  plot  it 

31  maxcols  =  max(J); 

32  maxrows  =  max(J'); 

[colmax,  xposition]  =  max (maxcols ) ; 

63 


33 


34 

35 

36 

37 

38 

39 


[ rowmax, 
colmax  = 
rowmax  = 


yposition]  =  max (maxrows ) 
colmax  -  noisefloor 
rowmax  -  noisefloor 


hold; 

plot (xposition,  yposition,  ' c+ ' ) 
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